Dynamics of coupled vortices in a pair of ferromagnetic disks 
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We here experimentally demonstrate that gyration modes of coupled vortices can be resonantly 
excited primarily by the ac current in a pair of ferromagnetic disks with variable separation. The 
sole gyration mode clearly splits into higher and lower frequency modes via dipolar interaction, 
where the main mode splitting is due to a chirality sensitive phase difference in gyrations of the 
coupled vortices, whereas the magnitude of the splitting is determined by their polarity configuration. 
These experimental results show that the coupled pair of vortices behaves similar to a diatomic 
molecule with bonding and anti-bonding states, implying a possibility for designing the magnonic 
band structure in a chain or an array of magnetic vortex oscillators. 
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Magnetic vortex structure is one of the fundamen- 
tal spin structures observed in submicron-sized ferromag- 
netic elements. It is well characterized by two degrees of 
freedom, one is 'chirality' (c = ±1), direction of the in- 
plane curling magnetization along the disk circumstance, 
and the other is 'polarity' (p = ±1), direction of the out- 
of-plane core magnetization. Particularly in the case of 
the disk, only the core region whose typical size is ~ 10 
nm generates the stray field, the static interaction be- 
tween vortices in an array is thus negligibly small in the 
ground state. However in the low frequency excitation 
state called the gyration (or translational) mode [H, 
the surface magnetic charges appear with the core mo- 
tion, that brings about the dynamic dipolar interaction 
between vortices. 

Coupled pair of magnetic vortices can be considered as 
a vortex molecule bound via dipolar interaction, which 
is a mimic of diatomic molecule with the van der Waals 
bonding [5j. The bonding or anti-bonding state respec- 
tively corresponds to in-phase or out-of-phase gyration of 
coupled vortices whose detailed energy levels are decided 
by combination of chiralities and polarities. This dy- 
namic coupling is also effective in a two-dimensional ar- 
ray system 0, Q , and thus allows us to design the density 
of states of the eigenfrequencies, the so-called "magnonic 
band structure", by arranging the core polarizations in 
a two-dimensional array. At the moment, there are few 
experimental reports for the coupled vortices via direct 
exchange interaction [8, 9] and also dipolar interaction in 
physically separated vortices However the prob- 

lem is still wide open in terms of experimental determi- 
nation of the detailed condition for the mode splitting 
and the magnitude of the dipolar coupling in interacting 
vortices. 

Herein we demonstrate the experimental evidence of 
the resonant excitation of a magnetostatically coupled 
pair of vortices as a clear mode splitting of spectra. 



The partial excitation using ac current causes the en- 
ergy transfer via dipolar interaction between two vor- 
tices, which results in a collective excitation of coupled 
gyrations. The observed mode splitting is reproduced 
by both micromagnetic simulation and analytical calcu- 
lation. The different combinations of core polarities and 
phase difference in gyrations give rise to four distinct 
modes. The magnitude of the vortex coupling is also 
discussed. 

Figure 1(a) shows a scanning electron microscope 
(SEM) image of the sample together with the measure- 
ment circuit. Each disk has the same dimension, 500 nm 
in disk radius r and 50 nm in thickness. Two neighboring 
Permalloy (FeigNigi; Py) disks with variable edge to edge 
separation d from 75 nm to 250 nm and electrical leads 
are fabricated on thermally oxidized Si (100) by means of 
electron beam lithography combined with electron beam 
evaporation techniques. The polarities of vortices were 
confirmed by means of magnetic force microscopy prior 
to all the electrical measurements. Hereafter, the num- 
ber 1 or 2 respectively represents the vortex confined in 
the left or that in the right disk shown in the SEM image 
of Fig. 1. One of the paired vortices is excited by a radio 
frequency ac current I ac T3-[l6l. [Til . fl9| and the result- 
ing dc voltage Vd c through a bias tee is synchronously 
detected by the same electrical contact probes [HI, |20| . 
In order to set the configuration of polarities P1P2 = 1 
or P1P2 = — 1, pi of the excited vortex is switched by 
applying high I ac of about 20 mA (3.5 x 10 11 A/m 2 ) at 
the resonant frequency [21] . 

Figure 1(b) shows the measured Vdc I Iac as a function 
of the frequency of I ac for the Py disk pair with d = 
75 nm. The reference spectrum for single disk shown by 
green symbols exhibits a sole dip at 352 MHz, which cor- 
responds to the resonant frequency of the vortex core gy- 
ration. When I ac is applied to one of the two neighboring 
Py disks, clear mode splitting takes place as can be seen 
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FIG. 1: (color on line), (a) Schematic diagram of the measure- 
ment circuit and an SEM image of the sample. Two Copper 
electrodes are attached to one of the disks in the Permalloy 
disk pair and the enveloped core is excited by a radio fre- 
quency current. Dynamics of the cores can be detected as dc 
voltages through the spin-torque diode effect utilized by a re- 
sistance oscillation associated with the core gyration. A lock- 
in technique is adopted at room temperature, (b) Frequency 
dependence of the normalized dc voltage Vdc/Iac measured 
for an isolated disk (green triangles) and for the paired disks 
with different polarities; black squares for pip2 — 1, and red 
circles for p\p2 = — 1. The ac current amplitudes I ac used for 
the measurements are I ac — 3.8 mA and 7.1 mA for the single 
disk, and I ac — 6.3 mA for the paired disks with the edge to 
edge distance d = 75 nm. Solid curves in each spectrum rep- 
resents the best fit to the data points using Eq. (1), thereby 
the dipolar coupling is evaluated. 



in the black and red spectra in the figure. Important to 
note is that the gyration mode of the single vortex with 
I ac = 7.1 mA only shows lower frequency shift and asym- 
metric broadening due to nonlinear effects [22| as plotted 
by open symbols. Therefore, the observed mode splitting 
is primarily due to the effect of the dipolar coupling be- 
tween the disks mediated by magnetic side charges. The 
magnitude of the mode splitting for P1P2 = — 1 is slightly 
enhanced by several MHz compared to that for P1P2 = 1. 

To gain insights into the dynamics of magnetostatically 
coupled vortices, micromagnetic simulations based on the 
Landau-Lifshitz- Gilbert equation [23[ were performed on 
pair of Py disks with identical physical dimensions. Typ- 
ical material parameters for Py are used : the saturation 
magnetization M s = 1 T, the exchange stiffness constant 
A = 1.05 x 10 -11 J/m, the spin polarization P = 0.4 and 
the damping coefficient a = 0.01. The disk is divided 
into rectangular prisms of 5 x 5 x 50 nm 3 for the simu- 
lation. A uniform I ac of 1.2 x 10 11 A/m 2 is applied only 
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FIG. 2: (color on line), (a) Simulated time evolutions of vor- 
tex cores at resonance frequencies for (piP2, c\Oi) — 
under ac currents (365 MHz and 390 MHz). Blue solid lines 
show the motions of the current-excited core in the left disk 
and red lines correspond to those of the indirectly excited 
core in the right disk, (b) Dispersion relations of amplitude of 
steady gyrations. Values of by max show the radii of steady gy- 
rations (50 nanoseconds after beginning of the current flow). 
Black symbols correspond to the parallel polarities p\p2 — 1 
and red ones to anti-parallel polarities p\p2 = —1 for oppo- 
site chiralities C1C2 = — 1. Results of same chiralities c\C2 — 1 
are plotted by open symbols. Simulation results for a single 
vortex are also presented by green symbols for comparison. 

to the left disk. 

After several nanoseconds from the start of excitation, 
the core gyration settles in an almost circular orbit and 
its amplitude is strongly enhanced at the resonance fre- 
quency [3, [Hj. This induces the core gyration in the 
neighboring disk, which also settles in the steady cir- 
cular orbit, and the collective gyration of the two vor- 
tices becomes fully synchronized and the eigenfrequen- 
cies of these modes appear in the spectra as character- 
istic resonance frequencies. Figure 2(a) represents the 
time evolution of the core deviations Sy at lower and 
higher resonance frequencies with respect to the single 
vortex for parallel polarities P1P2 = 1 and the same chi- 
ralities c\C2 = 1. At the lower frequency (365 MHz) 
both left and right cores rotate almost in-phase, whereas 
at the high frequency (390 MHz) the phase of the right 
core is retarded by approximately half a period. This is 
very much in analogy with covalent bonding in diatomic 
molecules or other forms of coupled oscillators. In the 
case of in-phase excitation, magnetic charges at side sur- 
faces form magnetic dipoles resulting in attractive force 
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between two cores, which corresponds to a bonding or- 
bital. On the other hand, the side charges of the disks 
repel each other for the out-of-phase excitation, having 
an analogy with an anti-bonding orbital. 

Figure 2(b) shows the maximum deviation of the core 
from the center (Sy max ) of the steady orbital for the cur- 
rent excited vortex in the left disk as a function of the 
ac frequency for the sample with d = 75 nm. The single 
vortex has a clear resonance peak at around 380 MHz 
and exhibits a small discrepancy with the experimental 
result (Fig. 1(b)) caused by self reduced magnetization 
(2lj . For the coupled vortices, two clear resonance peaks 
are observed on both higher and lower frequencies rel- 
ative to the sole peak for the single vortex. The mag- 
nitude of mode splitting for P1P2 = — 1 is larger than 
that for P1P2 = 1, as experimentally observed in Fig. 1 
(b). At both low and high frequencies, the in-phase and 
out-of-phase modes are degenerated with respect to the 
chirality. Figure 3 summarizes all four resonance modes 
from A to D characterized by rotational directions and 
the phase difference in gyrations. The lower frequency 
modes (A and C) stabilize with the help of attractive in- 
teraction due to the magnetic charges appearing along 
the disk circumferences. On the other hand, the higher 
frequency modes (B and D) stabilize with the help of re- 
pulsive interaction due to the opposite magnetic charges. 
Therefore the appearance of either in-phase or out-of- 
phase gyration in Fig. 2 (a) depends only on the polarity 
configuration piP2- 
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FIG. 3: (color on line). Schematic diagrams of four different 
resonance modes. Each mode is identified with the sign of 
polarities pi£>2, chirarities c±C2 and the phase difference of 
the gyrations. 

In both experimental and simulated results, the split- 
ting amplitude tends to be enhanced with the decrease in 
the normalized separation distance d n = d/r, as shown 
in Fig. 4. To check the effect of the current induced Oer- 
sted field on the excitation of the gyration of the right Py 
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FIG. 4: (color on line). Eigenfrequencies of the coupled gy- 
ration modes as a function of the separation distances. The 
simulation and experimental results are shown in (a) and (b) , 
respectively. The separation distance is given by a dimen- 
sionless value d n = d/R. (c) Estimated values of coupling 
strength rj x from both experiment and simulation with differ- 
ent separation distances. A d~ 6 curve from the rigid vortex 
model is also plotted as a dotted line. 



disk, the simulation was performed by replacing the left 
Py disk with an electrode where the current flows with 
the identical condition. The excited gyration amplitude 
due to the Oersted field is much smaller than that for 
the coupled vortices, implying that the magnetic dipolar 
interaction is the dominant factor for the indirect excita- 
tion of the right Py disk. 

We here discuss the experimental results analytically 
on the basis of Thiele's equation [UIHI, which describes 
the gyrations of vortex cores Ti = {X^Yi) with i = 1,2. 
It should be noted that our calculations take into ac- 
count only the adiabatic spin transfer torque as the ex- 
citation force of the vortex 1 in the left disk to simplify 
the analyses. The non-adiabatic and the Oersted field 
terms for the excitation are neglected. The effect of the 
non-adiabatic term fundamentally appears as the phase 
and the rotational radius shifts of the gyroscopic motion 
of the vortex core. One would point out the current in- 
duced Oersted field as the other origin of core excitation. 
Both the terms would affect the dynamics of the vortex 
1 however, may not influence the coupled modes 



identified in Fig. 3. Herein, the detected dc voltage can 



be evaluated by the following equation [20| ; 

V dc = 7 ac C/2(Re[X 1 ] cos^i + Re[Yi] sinJi), (1) 

where C is a constant and 5\ is the phase difference be- 
tween core and resistance oscillations in the left disk. The 
model [5] assumes the magnetostatic interaction energy 
term in Thiele's equation as 

U int = c lC2 /R 2 (VxXiX2 - n v YiY 2 ) + 0(\r/R 3 \). (2) 

According to the rigid vortex model [24| , the values of rj x 
and rj y are decided only by the shape of ferromagnetic 
element and the separation distance, independent from 
the excitation amplitude. As a simple solution for 
these problems, r] x and rj y should be treated as phe- 
nomenological fitting parameters including the influence 
of volume charges under the assumption of complete 
linear response with the fixed parabolic potential 
U = l/2k x X? + l/2k y Y? (k x ,k y = const.) and the 
interaction energy Ui nt . The best fitted curves to the 
experimental data using the above equation are shown 
in Fig. 1(b) by solid curves. The values of k x and k y 
are decided by a result of single vortex with I ac = 3.8 
mA as (k x , k y ) = (4.5 x 1(T 4 J/m 2 , 6.4 x 10~ 4 J/m 2 ). 
These results agree well with the previous work [20|, 
implying adequacy of Thiele's equation for the present 
analyses. The obtained values from experiments with 75 
nm separation are (rj x , rj y ) = (1.4 x 10 -17 J, 4.9 x 10 -22 
J) and those from simulations performed under identical 
conditions are (rj x , rj y ) = (1.2 x 10" 17 J, 8.1 x 10" 18 J). 
While the value of rj x agrees well each other, a large 
discrepancy of rj y is observed. This is possibly due to 
anisotropy along ^/-direction in the experiment due to 
the presence of attached Cu electrodes, which causes 
deviation of stray field from the simulation, resulting in 
closer resonance frequencies for coupled vortices with 
P1P2 = ±1 (red and black symbols) in Fig. 4(b). The 
separation dependence of rj x is shown in Fig. 4 (c). It 
is clear that the coupling strength increases@with@the 
decrease in@the@separation@distance, which 
can@be@clearly@reproduced by the two-dimensional 
micromagnetic simulation. It should be noted that the 
d dependence of rj x does not follow the d~ 6 dependence 
expected from a rigid vortex model {U @|. A strong 
magnetostatic coupling modifies the trajectories of the 
gyrations and the magnetization configuration near the 
edge region thus causes a deviation from the model, 
where a circular magnetization around the core is 
assumed. 

In summary, we have experimentally demonstrated the 
resonant excitation of coupled gyration modes in paired 
vortices by means of local excitation by an ac current 
passing through one of the disks in the pair. Excited 
coupled modes are identified by rotational directions and 
a phase difference as four different eigen- modes. The 



unique property in this system give us a guiding prin- 
ciple for designing the magnonic crystal in further ex- 
panded systems such as one-dimensional chains and two- 
dimensional arrays and a candidate for novel tunable os- 
cillators using vortices 0, 0, [27} . 
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